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Abstract 

We present the first study of the vector (Wess-Zumino) current in r~ — > K~ir~ K + u T decay using 
data collected with the CLEO III detector at the Cornell Electron Storage Ring. We determine the 
quantitative contributions to the decay width from the vector and axial vector currents. Within 
the framework of a model by Kiihn and Mirkes, we identify the quantitative contributions to the 
total decay rate from the intermediate states tvir, p^n and K*K. 
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Hadronic r decays provide a powerful tool to study low energy strong-interaction physics. 
The decay r~ — > K~n~ K + v T can proceed via both the vector and axial vector currents 
the vector current proceeds via the Wess-Zumino mechanism 0]. This model-independent 
mechanism plays an important role in hadron dynamics. It is applicable to r decays 0, |j| 
with three or more hadrons in the final state and violates the rule that the vector and 
axial vector currents produce an even and an odd number of pseudoscalars, respectively. 
Further interest in this decay is motivated by the fact that more precise knowledge of the 
KitK mass spectrum near the r mass is essential to give more stringent constraints on the 
r neutrino mass p|. In this Letter, we present the first study of the vector (Wess-Zumino) 
current, as well as the axial vector current, in the decay r~ — > K~7T~ K + i/ T , and determine 
their quantitative contributions to the decay width. Charge conjugate decays are implied 
throughout this Letter. 

In the Standard Model, the general form for the Cabibbo-allowed semileptonic r decay 
matrix element can be expressed as M. = ■%V u du(q u )'~f fl (l — , jF J )u(q T ) J M , where G is the Fermi 
coupling constant, V u d is a Cabibbo-Kobayashi-Maskawa matrix element; q v and q T are the 
four- momenta of the r neutrino and the r lepton, respectively; = (hadronsl^^ — v4 M |0) 
is the hadronic current, and V 1 and A^ 1 are the vector and axial vector quark currents. 
The most general ansatz for the hadronic current of three hadrons is characterized by four 
form-factors |j| in terms of the hadrons' momenta q^ (here, i — 1 for K~ , i = 2 for n" and 
i = 3 for K+), 

J" = (q^-q$-Q»^^y^s^) 

+ ^-^-Q" 9 ^^)f 2 (s 1 ,s 2 ,Q 2 ) 
+te^q la q 2 ^F 3 (s u s 2 ,Q 2 ) 

+Q fl F i (s 1 ,s 2 ,Q 2 ), (1) 

where = + q 2 + q% . All four form-factors F\ through F± are functions of Q 2 = 
(<?i + <?2 + q 3 ) 2 , s 1 = (q 2 + g 3 ) 2 , s 2 = (qi + q 3 ) 2 and s 3 = (qi + q 2 ) 2 . The terms proportional to 
Fi and F 2 originate from the axial vector current (J p = 1 + ). The vector current (J p = 1~) 
originating from the Wess-Zumino mechanism gives rise to the term proportional to F 3 . F 4 
is due to the spin-zero scalar current (J p = + ) that is expected to be small [3,|f| and in this 
analysis it will be set to zero, although we will consider its contribution to the systematic 
uncertainties later. 

Tau decay to three mesons is conveniently analyzed in the hadronic rest frame, where qi+ 
q 2 + q 3 = 0. The angle /3 @ is defined as the angle between the direction of the hadronic 
system in the laboratory frame and the normal to the plane defined by the momenta of 
particles 1 and 2, 

cos/3 = — Q • norm(qi x q 2 ). (2) 

The angle 9 0, defined to be the angle between the flight direction of the r lepton in the 
laboratory frame and the direction of the hadronic system as seen in the r rest frame, is 
related to the energy Eh of the hadronic system in the laboratory frame by 

2xm 2 T — m 2 — Q 2 
cosO = , =, (3) 

{ml - Q 2 )Jl - Am 2 r /s 
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where x = 2Eh/\fs and s is the square of the center of mass energy. The angle if) between 
the flight direction of the r lepton and that of the laboratory frame as seen from the hadronic 
rest frame, — Q, is given by 

. x(ml + Q 2 ) - 2Q 2 

cosw = j ■ (4) 

(m 2 -Q 2 ) Jx 2 -AQ 2 /s 



The differential decay width for the decay r — > K tt K + v t with the scalar contribution 
neglected, after integrating over the unobserved neutrino direction, is given by 



dT(r -> KitKv T ) G 2 2 1 (m 2 T -Q 

\Vud 



2\2 



dQ 2 d Sl ds 2 12m T 1 uai (4tt) 5 Q 4 

+ + (5) 

Here, the structure functions Wa and Wb can be expressed in terms of the form factors as: 

W A = (x 2 + x 2 )|F 1 | 2 + (x 2 + x 2 )|F 2 | 2 

+2(x 1 x 2 -x 2 )Re(F 1 F 2 *), 
W B = x 2 |F 3 | 2 , (6) 

where Xi jl| are known functions of Q 2 , Si, s 2 and S3. Thus and W^b govern the rate 
and the distributions of Q 2 , S\, s 2 and S3. There is no interference between the axial vector 
current contributions (F 1 ,F 2 ) and the vector current contribution (F 3 ) to the decay width. 

A parameterization of the form-factors is given in |6j], based on chiral perturbation theory 
at low momenta and meson resonance dominance at higher momenta. It has been imple- 
mented in the Monte Carlo program koralb |7| as: 

Fi = _V2 BWai(g2 ,BW p (s 2 ) + ^BW,(s 2 ) 



V2 „ , o2l 



p 



F 2 = -fj- ■ Rf -BW ai (Q 2 ) -EWk^si), 

<->J-K 

1 r— BW u (s 2 ) + aBWjf.(si) 



2\/27r 2 /3 v 1 + a 

BW P (Q 2 ) + ABW P '(Q 2 ) + (5BW P »(Q 2 ) 

where / ff = 0.0933 GeV and f3 p = -0.145 0,0- BW x (s) is the two- (dependent on s x or 
s 2 ) or three- (dependent on Q 2 ) particle Breit-Wigner propagator with an energy dependent 
width. The CLEO result on the a± parameterization derived from t~ — > 7r _ 7r°7r°z/ r decay 0] 
is used. The parameterizations of two-particle Breit-Wigner functions take the forms used 
in . The parameters used for the three-particle Breit-Wigner functions are taken from the 
Particle Data Tables @: m p = 0.770 GeV, T p = 0.151 GeV, m p , = 1.465 GeV, T p > = 0.310 
GeV, m p n = 1.700 GeV and I> = 0.240 GeV. 

The form-factors F\ and F 2 correspond to the axial vector current processes a\ — ► p^ir~ 
with p('> — ► K~K + and ai — > K*K~ with if* — > if + 7r~, respectively. The form-factor F3 
represents the vector current processes decays to K*K~ and u;7r~ with K* — > if + 7r~ and 
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uj — > K K + . G parity conservation forbids the vector current process W — > pV'"} — > p^n 
with pW K~K + . Instead, in @ the process W — ► p^'"-* — > cj7t is considered, although 
the subsequent decay uo — > K~K + proceeds only through the high- mass tail of its Breit 
Wigner propagator. We make use of this model in this analysis, although we will be unable 
to distinguish whether the decay proceeds via the uj meson or via some other resonance with 
pole mass outside the kinematic limits for the KK mass. 

The default version of koralb |?J does not model the decay r~ — > K~n~ K + v T well, 
as indicated by the discrepancies in describing the invariant mass distributions [10j . To 
improve the agreement, we introduce two more parameters Rb and Rp, describing the 
relative strengths of Wb and Wa and of F 2 and Fx in Eqs. |S] and El These two parameters 
are explicitly expressed in our generalization of the form-factors in Eq. [71 The five model 
parameters: a, A, 5, Rb and Rp are determined in this analysis by a simultaneous fit to 
the invariant mass spectra of K~tc~K + , K + ti~ and K~K + . These five parameters allow for 
a fit that properly describes all the mass spectra. In general, these parameters should be 
complex, allowing additional interferences. However, due to limited statistics, we have set 
them real and taken into account only the inherent phase of the Breit- Wigner functions in 
Eq.UJ 

The data sample used for this analysis corresponds to an integrated luminosity of 7.77 fb _1 
taken on or near the T(45). It contains about 7.09xl0 6 tau-pairs. The data were collected 
with the CLEO III detector located at the e + e~ Cornell Electron Storage Ring. The 
CLEO III detector configuration features a four-layer silicon strip vertex detector, a wire 
drift chamber, and, most importantly for this analysis, a ring imaging Cherenkov (RICH) 
particle identification system. A detailed description of the RICH performance can be found 



in 



0. 



To select the three-prong r decays to the K~n~ K + u T final state, we use one-prong decays 
of the other r to e + h>v, p + vv, p + v and -k + {K + )v to tag tau-pair events. Pion and kaon 
identification is obtained by a combination of the RICH information with dE/dx measured 
in the drift chamber. We use the event selection criteria listed in and find 2255 candidate 
signal events in our data sample. The branching fraction obtained using this data sample 
for the decay t~ — ► K~ix~ ~K + u T is consistent with the previously published result |13j |. 

Background contributions and efficiencies are obtained using Monte Carlo events from 
the koralb [7] and jetset w]| generators for the tau-pair production and qq — > hadron 



processes, respectively. These events are then processed by the GEANT-based [15J] CLEO 
detector simulation and pattern reconstruction. In our sample, 256 ± 16 ± 46 events are 
attributed to the backgrounds from the continuum qq production and from tau decays to 
other channels. The dominant backgrounds are due to pion mis-identification or missing 7T° 
from the tau cross-feed decays t~ — ► K~7t + tt~v t , 7r~7r + 7r~z/ r , and K~tt~ K + tt°v t , or from 
t~ — > K~7i~ K + v T decay with a mis- identified tag side. Each of these sources contributes 
about 20% of the total backgrounds. The events with the three-prong side identified as 
K~7T~ K + u T and the tag side mis-identified are analysed in the same way as the signals. 

We use the unbinned extended maximum likelihood method to extract the parameters 
of the hadronic structure of the decay described above. The likelihood function, C, has the 
form 

e~ m m n n 1 

C = j— Y[-\N s V^(a,X,5,R B ,R F ,Q 2 ,sx,s 2 ) 

+ N B V?( aj A, 5, R B , R f , Q\ Sl , s 2 )} . (8) 
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The number of signal and background events are denoted by N s and N B with n = N s + N B . 
m is the number of the expected events. Vf and Vf represent the normalized probability 
density functions (PDF) for event % to be either a signal or a background event. The 
signal PDF is formed from the product of one- dimensional PDFs obtained from integrations 
over the differential decay width described by Eq. |5J convoluted with the mass-dependent 
efficiencies. The background PDFs are obtained from Monte Carlo simulations. The fitting 
procedure has been tested with Monte Carlo samples to verify its performance. 
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FIG. 1: Projections of the fits (full lines) onto the a) K-ir~K + , b) K+tt' and c) K~K + mass axes 
from the data (dots with error bars) superimposed with the contributions from the vector (Wess- 
Zumino) current (dashed) and the axial vector current (dotted) and all backgrounds (dot-dashed) . 
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FIG. 2: Angular distributions of a) cos/3, b) cos# and c) cos^ for the data with statistical er- 
rors only. Solid lines represent our model of the vector (Wess-Zumino) and axial-vector signal 
contributions and all backgrounds, using the parameters from the full fit to the data. Dashed 
lines represent our model of signal and backgrounds with the vector (Wess-Zumino) contribution 
excluded. Dotted lines show all backgrounds. 



The results of the unbinned extended maximum likelihood fits are projected onto the 
K~-k~K + , K + ir~ and K~ K + mass axes with different contributions superimposed as illus- 
trated in Fig. We observe that about half of the r decay to K~n~ K + v T proceeds via the 
vector (Wess-Zumino) current. The fit parameters are a = 0.471 ±0.060, A = — 0.314±0.073, 
5 = 0.101 ± 0.020, R B = 3.23 ± 0.26 and R F = 0.98 ± 0.15 (statistical errors only). The 
angular distributions of cos/5, cos 9 and cos^ for the data and for the two Monte Carlo 
samples generated using the parameters from the fits are shown in Fig. El The comparison 
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shows that the Monte Carlo samples with and without the contribution from the vector 
(Wess-Zumino) current both describe the data well within the current statistics; the angular 
observables therefore have little sensitivity to the presence of the vector (Wess-Zumino) cur- 
rent in our sample. Our ability to distinguish the axial vector contribution from the vector 
(Wess-Zumino) contribution comes from the invariant mass distributions, and thus depends 
on the correct modeling of the associated form-factors, as constrained by our data. 

We have considered the following sources of systematic errors: the systematic errors for 
the five model parameters and for the fractions of the vector and axial vector current con- 
tributions to the decay rate (including those of the intermediate states) come from the mass 
dependence of the detection efficiency, the uncertainties in the modeling of the mass resolu- 
tions, the uncertainties in the parameterizations of the resonances and from the neglected 
contribution of the scalar current. To estimate the effects of the detection efficiencies, we 
introduce an artificial linear mass dependence of the efficiencies and vary the corresponding 
slopes by ±10%. We estimate the associated systematic error on the relative strength of 
the vector and axial vector currents as ARb/Rb — 10%. The corresponding systematic 
errors on all other fit parameters and the fractions of the vector and axial vector current 
contributions to the decay rate are less than 4% of their values. The typical mass resolutions 
of the hadronic final state in the decay r~ — > K~ir~K + v T are about 2-3 MeV, i.e., much 
smaller than the widths of the resonances considered. The uncertainty in the modeling of 
the detector resolutions contributes a <2% systematic error. The uncertainty introduced 
by the neglected scalar current is estimated at the model predicted level p, 0] and found 
to be less than 1%. In addition, when we vary the normalization and the shapes of the 
background PDFs by ±ler, the fit results vary between 1% and 12%. When the parameter 
j3 p in Eq. [7| varied by ±100%, the results change 37% for the parameter R F and less than 
3% for all other fit parameters and the contributions of the vector and axial vector currents. 
The dominant uncertainty in the fit results arises from the errors of the parameterizations 
of the poorly measured p' and p" resonances 0. The corresponding contribution to the 
overall systematic uncertainties for the parameters a, A, 5, Rb, and Rf is 5%, 23%, 153%, 
56%, 1%, respectively. Similarly it is 8%, 10%, 11%, 12%, 9% and 27% for the fractions 

TZwz, T^av K , T^av 7 "' ^wz^ aric ^ ^wz to be presented. We add all estimates of the 

systematic errors in quadrature and obtain the overall systematic uncertainties for the five 
parameters a (7%), A (25%), 5 (154%), R B (59%), and R F (37%) and for the fractions U wz 

(9%), TZ AY (H%), Tl% K (11%), TZ P ^ (16%), n«* z K (10%) and TZ^ Z (29%). There is a 
large uncertainty in the contribution of the p" resonance to the vector current. This has 
little effect on the total vector contribution because it is strongly suppressed by phase space, 
and its contribution is easily absorbed into the contribution from the pf resonance. 

In summary, we have presented the first study of the vector (Wess-Zumino) current in 
the decay r~ — > K~tt~ K + u t and determined its contribution, as well as that from the axial 
vector current, to be: 

TZ WZ = ^ = (55.7 ±8.4 ±4.9)%, 

1 tot 

TZ AY = ^ = (44.3 ± 8.4 ± 4.9)%. (9) 

Itot 

The errors are statistical and systematic, respectively. 

The structure of the decay r~ — > K~tt~ K + v t can be modeled using the parameterization 
given by Eqs. |3HZ1 with the parameters determined to be: a = 0.471 ± 0.060 ± 0.034, 
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A = -0.314 ± 0.073 ± 0.080, 5 = 0.101 ± 0.020 ± 0.156, R B = 3.23 ± 0.26 ± 1.90, and 
R F = 0.98 ±0.15 ±0.36. 

In the context of the model by Kiihn and Mirkes (if, where the axial vector current 
proceeds via a\ — > p^n~ and a\ — > K*K~ , and the vector (Wess-Zumino) current via p(''"' — > 
K*K~ and ujtc~ , we can also separate the individual contributions from the intermediate 
states ujix~ , p^ix^ and K*K~ to the total decay width as 

= = (3.4 ± 0.9 ± 1.0)%, 



p(')7r 



p (/) 7r _ r^y 



Ufa* = = (2.5 ± 0.8 ± 0.4)%, 

1 tot 

K$l K = = (60.8 ± 8.5 ± 6.0)%, 

Ttot 

n x*K = Ejv^ = (46>8 ± 8 4 ± 5>2)%> (10) 

1 tot 

These fractions do not add up to 100% due to the interference of the amplitudes of the 
intermediate states. The decay is dominated by the vector (Wess-Zumino) and axial vector 
current processes via the intermediate state K*K. The ratio R%j K , together with Bir~ — > 
K-n-K + v T ) H, yields B(at K*K) = (2.2 ± 0.5)% which is compatible with the result 
B(ai — > K*K) = (3.3 ± 0.5)% extracted from the decay r~ — > 7r~7r°7r t V T j8j with only one 
threshold (K*K) considered. The apparent shortfall of the 7r~7r°7r° mass spectrum in the 
previous analysis is attributed to the neglected contribution from the vector (Wess-Zumino) 
current. The axial vector current contribution determined in this Letter is much lower than 
the result from ALEPH ^| based on CVC and limited e + e~ data. As noted above, we 
can not distinguish the p(' , "> — > tun~ with uj — > K~K + decay from other models, such as 
other vector meson resonances or even a simple constant. Setting the ujh~ contribution to 

zero, we obtain Tl^ K) = (50.8 ± 7.7)%, TZ AV = (49.2 ± 7.7)%, TZ P ^ = (4.5 ± 1.4)%, and 
TZ^y K = (51.6 ± 7.8)%. However, the fit favors the presence of an additional resonance 
component, which can be associated, within the context of the model in 0, with utt~ , over 
its absence by nearly Qa (statistical error only). The model parameters presented here can 
reduce the model dependent uncertainties on the tau neutrino (z/ r ) mass measurements using 
the decay r~ — > K~ix~ K + v T . 
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luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, and the Texas Advanced 
Research Program. 
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